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Abstract 

Wc show that at asymptotically high densities the "color-flavor-locked + neutral kaon condensate" phase of 
QCD develops a charged kaon condensate through the Coleman- Weinberg mechanism. At densities achiev- 
able in neutron stars a charged kaon condensate forms only for some (natural) values of the low energy 
constants describing the low-lying excitations of the ground state. 
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It was realized a long time ago |T], ^, ^ that the attraction between two quarks close to the Fermi 
surface in high density strongly interacting matter leads to the formation of Cooper pairs of quarks 
and the spontaneous breaking of color symmetry. This phenomenon was more recently studied using 
Nambu-Jona-Lasinio models and renormalization group methods Q ^, ^, instanton models ^ 
and, at asymptotically high densities where a weak coupling expansion is valid, perturbative QCD 
1^, The main lessons learned in these studies were that i) the gap can be large, up to 100 

MeV, ii) in the case of three-flavors of quarks with the same mass the ground state is the so-called 
color-flavor-locked (CFL) state , and chiral symmetry remains broken at arbitrarily high densities 
and iii) the low-lying excitations carry the quantum numbers of the pseudo-scalar octet familiar 
from the zero density case (plus two other scalars related to the spontaneous breaking of baryon 
number and axial charge) and iv) a form of electromagnetism survives: a combination of the photon 
and the eight gluon is not "Higgsed" and remains massless in the CFL phase. The equation of state 
is left nearly unchanged by the pairing. A number of electromagnetic and transport properties of 
quark matter are, however, sensitively dependent on what phase the system finds itself and what 
the low energy excitation are. This dependence provides a unique opportunity to study quark 
matter in the interior of neutron stars (or rule out its existence). 

The fact that the low-lying excitation of the ground state are very similar to the low-lying 
excitations of the vacuum (pions, kaons,...) allows us to studied small perturbations around this 
ground state with the techniques of chiral perturbation theory. This brings about two advan- 
tages. At asymptotic high densities where perturbative QCD is valid it organizes perturbative 
calculations that would be very complicated otherwise. More importantly, it provides a method 
to systematically expand around the CFL phase in inverse powers of the density and/or the gap, 
which is particularly useful if the use of perturbative QCD is not legitimate. One of the impor- 
tant perturbations around the CFL phase is the presence of realistic quark masses. In the case 
of free quarks it is easy to determine the response of the system to these masses. As the mass 
of the strange quark increases towards its realistic value its density decreases in such a way that 
its Fermi energy equals the Fermi energy of the up quark (plus the Fermi energy of the electrons) 
and the weak decay s ^ u + + I'e becomes forbidden. Charge neutrality is guaranteed by the 
presence of electrons (we assume that the neutrinos leave the system). The interacting case may 
be qualitatively different. If the interactions are such that quarks of different flavors are paired, the 
change of flavor caused by flavor changing decays would result in two unpaired quarks, what is not 
energetically favorable. As a result, the system is rigid against small enough flavor asymmetries. 
That is what happens in two-flavor QCD, where up and down quarks (of two of the three colors) 
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are paired. An asymmetry in mass or chemical potential between the flavors causes little change 
in the ground state if they are small enough, the small change coming entirely from the unpaired 



quarks of the third color |12]. One might think that the same effect occurs in the three-flavor, 
CFL phase, since there all quarks are paired with quarks of different flavors. The CFL phase 
has, however, a another way of responding to mass asymmetries that costs little energy but is not 
available in the free or two- flavor system : it can condense mesons carrying strangeness, that are 
particularly light This was demonstrated in on very general grounds and, in the case 

of weak coupling, through explicit computations of the response function to mass asymmetries. 
For realistic values of the quark masses and densities it was found that the is the meson that 
condenses and this ground state will be referred to from now on as the "CFL + K^" phase. For 
more general asymmetries, including asymmetries on chemical potential present before weak equi- 
librium is achieved and/or neutrinos leave the system, a rich phase diagrams results, with kaonic, 
pionic, neutral and charged condensates forming with different values of the parameters [19|. The 
pattern of symmetry breaking caused by the condensation (in the isospin limit) is the same one 
found in the standard electroweak model SUj{2) x i7y(l) — > Uq{1) (/ =isospin,y =hypercharge 
and Q = modified electric charge). Due to the lack of Lorentz symmetry, only two, and not three, 
Goldstone bosons are generated, one neutral and another charged |^, There are stable, 



superconducting topological vortices [22, 23 1 and almost stable non-topological domains walls in 
the "CFL + K^" phase [^]. In the CFL phase there is an equal number of quarks of the three 
flavors, and the system achieves electrical neutrality in the absence of any electrons, making it a 
perfect insulator |^5|. The presence of the condensate, being neutral, does not change this 
situation. A charged kaon condensate however would change quark matter from a perfect insulator 
to a (electrical) superconductor. It is a generic feature of charged massless scalars that the strong 
long wavelength fluctuations of the gauge field lead to condensation of the scalar field (Coleman- 
Weinberg mechanism (||2^)). In the CFL+K*^ phase there is one almost massless charged scalar 
field. Its mass comes from isospin breaking contributions coming from the quark mass difference 
and electromagnetic mass effects. In this paper we consider the competition between the isospin 
breaking mass terms and the fluctuations of the electromagnetic field in order to determine the 
fate of the charged kaons and of the possibility of a (electromagnetic) superconducting phase in 
quark matter. 

In the absence of quark masses, the symmetry breaking pattern generated by diquark conden- 
sation in the CFL phase is § SUc{3) x SUl{S) x SUr{3) x Ub{1) x Ua{1) ^ SUc+l+r{3) x ^2- 
The electromagnetic Uq{1) is a subgroup of the chiral group 5C/l(3) x SUr{3) and the there is a 
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surviving local Up^{l) "electromagnetism" in S'C/c+L+i?(3) that is a combination of the photon and 
one of the gluons. The axial Ua{^) is only an approximate symmetry of high density QCD due to 
the instanton suppression in the medium. This symmetry breaking pattern implies the existence 
of two singlet Goldstone bosons associated with the broken baryon number and axial symmetry, 
and an octet of pseudoscalars. The singlets will not play a role in our analysis and will be dropped 
from now on. At low (excitation) energies below the gap, QCD is equivalent to the most general 
theory of an octet of pseudoscalars and photons with the same symmetries of QCD. This theory 



has been extensively analyzed |27, 29, ^ recently. The leading terms of its lagrangian are 



e - 1 - aA^ ~ 

C = -E"^ - -B'^ + ^Tr[L>oS'^DoS - i^^VS+VS] + Tr[A^(S + S^" - 2)] 

2 2 4 Svr 

+ 6d/2A2Tr[St,Q][Q,S]+cd2/' (Tr[St,Q][Q,S])2 + .... (1) 

where ^ is the baryon chemical potential, A is the gap, DqTi = d^T, + ^[A^^,S] — ieAo[(5,S], 
i5s = VqS — ieA[Q,li\, e = eg\f?> j \/ ?>g^ + 46^ (gr =strong coupling constant, e the electron 
charge), a = e^/47r, }A and Q are the quark mass and charge matrix, }A = Al~^detA^ and 
T, = e . A few comments are in order here. The electromagnetic field in Eq. (Ill) are the 

rotated fields that remain massless in the CFL phase. The low energy constants /, v, a, b, c, e can, 
in principle, be determined from QCD. In practice, this can be done only in the asymptotic limit 
where perturbation theory is valid. At lower densities one can estimate their values by looking 
at their variation with the cutoff of the effective theory (see below). The dielectric constant e 



was computed in [27| where it was found that e = 1 + ^^^t- The magnetic permeability was 
argued to be unchanged from the vacuum value because the diquark condensate carries no spin. 
We will assume this to be true even outside the perturbative QCD regime. The values of /, v were 
also determined in perturbation theory to be f = l/\/3 and / = (21 — 81n(2))/367r^. After 
some controversy |3^, 32, 33] the value of a seems to have settled at a = 12/7r 33]. 



Finally, the value of the gap is estimated to be around 50 — 100 MeV in phenomenological models 



at ;U ~ SOOMeV and is given in perturbation theory by A = 512(2/ NF)^/27r'^/g^e J|, |T^, 

although large corrections from higher orders are expected ]|T^ , [T5| , [l6|] . The coefficient of the 
term quartic in M in Eq. (0) is not a free parameter because it is related by an approximate local 
symmetry of high density QCD to the kinetic term |13]. Terms violating this symmetry (like the 
other mass term in Eq. (||)) are suppressed by extra powers oil/ ^i. The electromagnetic coefficients 
b and c have not yet been computed in perturbation theory but we will estimate them below. The 
terms implied by the dots in Eq. (|l]) are further suppressed by powers of momenta or meson masses 



in units of the cutoff A ~ 2A or extra powers of a. 
For values of A satisfying 



Y_ J _^__u ajs ^^2 J s _ 

ms{ma-mu) 



the minimum of the potential is found at |13] 

/ 
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(3) 



with 



V2(l) am„AV 
/ TTf-^^ms - mdY(ms +rnd) 

describing a condensate (for reasonable values of the parameters cot is nearly one). The 

upper limit in Eq. (|2|) is the maximum value of A for condensation and the lower limit marks 
the onset of condensation. At very high /i the numerical value of the range in Eq. is fairly 
independent of the chemical potential ^i. Taking the perturbative QCD values of /, a, v and e, 
niu = 4 MeV and = 150 MeV this range is (2MeV)2 < A2 < (120MeV)2) and thus, most likely, 
the real world case will correspond to the "CFL+K'^" phase. Notice that in the CFL phase the 
kaons are the lightest mesons and so are the first to condense under perturbations. In the isospin 
limit neutral and charged kaons are degenerate. The reason for the neutral kaons to condense 
are their slightly smaller mass due to the — niu quark mass difference and electromagnetic 
corrections and the fact that the presence of a charged kaon condensate implies the presence of 
electrons to guarantee charge neutrality, which raises the energy of this state compared to the 
CFL+KO phase. 

This conclusion may be changed by the inclusion of photon loops and electromagnetic interaction 
terms that, by consistence, must be included together. Let us now compute the one photon loop 
contribution to the effective potential and leave the discussion of the conditions under which it is 
important for later. 

It is convenient to use a modified Landau gauge fixing procedure, that is, we add a term 
— l/2^((9o^^ + f^Vyl)^ to the lagrangian, taking the limit ^ — > 0. The price payed by having a 
complicated propagator that breaks Lorentz invariance is compensated by the fact that all zero 
external momentum one-loop diagrams involving a meson propagator vanish, since the photon 
propagator satisfies, in this gauge 



FIG. 1: Graphs giving rise to the one loop cfFective potential in the Landau gauge. Solid lines arc mesons, 
wiggly line are photons. 



P^V^'D^xip) = 0, V'^'' = diag(l, v', v\v'). 



(5) 



The one- loop effective potential is then given by the sum of diagrams shown in Fig.( 
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where we work in D = 4 + r/ dimensions and A is an arbitrary renormahzation scale. The ultraviolet 
divergences are absorbed in the terms proportional to h and c in Eq. (|l|), which suggests that the 
natural values for these (renormalized) constants at the cutoff scale A ~ 2A are 



b{\ ~ 2A) 



h ,2v'^ 1 , 
-(— + —. 



c(A=.2A) = ^(^ + 4)- 

16 y/e 



(7) 
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where b and c are numbers of order one^. In addition we expect on physical grounds that 6 < 0, 
what guarantees a positive electromagnetic contribution to the mass square for the mesons (the 
photon loop contribution vanishes in the Landau gauge). This estimate of the coefficient b agrees 
with the ones in |l^ (where no attempt was made to count factors of Att) and [^] (where no 
attempt was made to count factors of e or t;). Other electromagnetic terms not renormalized at one 
loop order are assumed to be suppressed. Also, the contribution from meson loops is proportional 
to (mK/f)' ~ (A^)4((m u + 77i^)7Tis/ /)^ and is strongly suppressed. Finally, the effective potential 
including the one photon loop correction becomes 

^-// = -^TVA^(S + St-2)--^TV[St,M2][_^2^q 
Svr 

- 6f^ + Al^Trpt,Q][Q,S] (8) 



In 



+ (^ + 4)^(Tr[St,Q][Q,S]) 



i;2e2/2Tr[Et,Q][Q,£] \ _ 3 
8A2 2 



Whether the electromagnetic corrections computed above can modify the position of the minimum 
of the potential depends on the hierarchy assumed for the scales fi, A, etc. 

It is probably instructive to compare the present situation with the simpler one of scalar QED. 



The effective potential there has the form (omitting numerical factors) [26| 

/,2\ 



y = mV + A0^ + (A2 + aV(^l- (9) 



1 , A 



In the massless case the minimum of (^ is at = M'^/a e~2+^. Assuming A ^ term 
proportional to can be disregarded and higher loop corrections are under control at cj) = cj). 
The presence of a finite mass term will not destroy this minimum if, at (j) = (p, it is smaller than 
the other terms. This condition translates into m?e~^~^'^ < aM"^. In our case, we have mass 
terms of the order aA^, (m^ — 'mu)msA.'^ / fJ.'^ , and A^. The role of the self-interaction 

is played by terms coming from the electromagnetic mass term (A ~ aA'^/ fi'^), the quark mass 
term (A ~ (jn^ — mu)msA'^ / fi"^) and the electromagnetic interaction term (A ~ a'^). Assuming 
a ~ A'^ / fi^, {nid — mu)ms/ (or larger) the condition for the survival of the non-trivial minimum 
is satisfied up to numerical factors. Those numerical factors (depending, among other things, on 
the low energy constants b and c) determine whether the Coleman- Weinberg mechanism occurs. 
Thus, let us analyze them more carefully in two separate situations. 



^ The values of h and c are gauge and renormalization prescription dependent. We refer here to their values in the 
Landau gauge ^ and in the modified minimal subtraction scheme with renormalization scale A ~ 2A 
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FIG. 2: £^ffective potential as a function of a; = t&TL{\K+\/\K°\) = 10 GeV, b = -1.3, c = 0.3, and leading 
pcrturbativc results for the remaining parameters). The dashed line shows the effective potential without 
the electromagnetic contribution. 

Asymptotic limit: To analyze the simultaneous condensation of neutral and charged kaons we 
use the parameterization 

K+ = (psinxe'^^. (10) 

Due to charge and hypercharge invariance the potential does not depend on the phases 6i and 62. 
Isospin breaking effects create a dependence on x. The effective potential (^ becomes 

Veff = {ma - rriu) ( " ^"^^^ I cos2x - — (— + f X 

+ + In ( + oi^JI^^^) (11) 



16 V8A2e^+^ / rn. 



s 



with X = 2sin2 ^(cos^ x + {1 + cos(-^^)) sin^ x) and we approximate cos(-^y^) ~ 0. 
Let us now consider the limit /i — > 00. For values of fi such that 

5m nis < (12) 

(but still satisfying the condition for K'^ condensation in Eq. (|^, that is violated only around 10^ 



GeV) the first term in Eq. (11) is smaller than the second one and can be disregarded. Numerically, 
condition ( p!2|) is satisfied for fi > 3GeV. Let us momentarily put aside the second the term in 
Eq. ( p!Tl) (electromagnetic mass). Minimizing in relation to x we find a solution 
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FIG. 3: £:ffective potential as a function of a; = tan(|i^+|/|i^O|) {^l = 500MeV, A = 50MeV, / = S/pew, & = 
—0.5, c = 2.5 (lower curve) and 1.5 (upper curve) and leading pcrturbativc results for the remaining param- 
eters). The dashed line shows the effective potential without the electromagnetic contribution. 

At this value of x, the electromagnetic mass term can be disregarded compared to the one we 
kept if 

- 46 e-(=+t) < 1. (14) 
For many, but not all, natural values of 6 and c this condition is satisfied and the solution in 



Eq. (13) can be trusted. Unfortunately, the asymptotic values of these parameters in the limit 



— > oo are not known (the computation of c involves the calculation of four loops diagrams) and 



we cannot determine whether (14). In Fig.(g) we show, as an example, the effective potential for 
a natural choice of parameter values and very high value of the chemical potential (/x = 10 Gev, 
rriu = 4 MeV, = 7 MeV, = 150 Mev, b = —1.3, c = 0.3). It shows the characteristic shape of 
a potential with a first order phase transition. 

"Realistic" densities: For fi < 3 Gev the quark mass terms are no longer negligible compared to 
the electromagnetic mass terms. In fact, for the densities that may be found in neutron star cores 
(/U ~ 500 MeV) it is the dominant mass term for the charged kaons and the one loop effects are 
too small to overcome it for most values of the parameters. However, at lower densities the values 
of the low energy constants are not so well determined since the perturbative results do not apply. 
Some choices for the values of these low energy constants that do not violate the expectations of 
dimensional analysis result in charged kaon condensation. As an example we show in Fig. (I) the 
effective potential for two choices of the parameters. 

In both of them the value of the decay constant / was changed from the value suggested from 



perturbation theory (/ O.Gfj. = 3 fpert) and used [b] = —0.5. The two sohd curves correspond 
to c = 2.5 and c = 1.5. This change in c is enough to transform the global minimum into a local 
minimum. This set of parameters were carefully chosen For most of the parameter space the quark 
mass term overwhelms the others and there is no charged kaon condensation. 

A better idea of the likelihood of charged kaon condensation at these densities can be perhaps 
obtained through the use of QCD models to estimate the unknown low energy constants in the 
density range inaccessible to perturbation theory. 

We have considered the possibility of charged kaon condensation and (electromagnetic) super- 
conductivity at high dense quark matter. At asymptotically high densities, where perturbative 
QCD applies and the question can be decided on first principles, a complicated computation of 
some low energy constants are necessary to settle the issue. We find however that for most natural 
values of these constants charged kaon condensation indeed occurs. At lower densities the situa- 
tion is the opposite. For most reasonable values of the low energy constants the quark mass effects 
overwhelm the electromagnetic effects and there is no condensation. 
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